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Abstract 


This note studies classical magnetohydrodynamic shock waves in an inviscid fluidic 
plasma that is assumed to be a perfect conductor of heat as well as of electricity. For 
this mathematically prototypical material, it identifies a critical manifold in parameter 
space, across which slow classical MHD shock waves undergo emergence of a complex 
conjugate pair of unstable transverse modes. In the reflectionally symmetric case of 
parallel shocks, this emergence happens at the spectral value A = A/|a;| = 0, and the 
critical manifold possesses a simple explicit algebraic representation. Results of refined 
numerical treatment show that for only almost parallel shocks the unstable mode pair 
emerges from two spectral values A = ±^ 7,7 > 0. 


1 The equations of ideal isothermal MHD 

We consider ideal MHD in two dimensional space, 

0 =pt + div(pV) 

0 ={pV)t + div(pV ®V + {p+ \\B\^)I -B®B) (1) 

0 =Bt + div{B ®V - V ® B). 

The dependent variables p > 0,p > 0, V G denote the flnid’s density, pressnre, and 
velocity. In addition to Q, the magnetic held i? G satishes 

divH = 0 . ( 2 ) 
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The fluid is assumed to be polytropic, p = RpT, and have a constant temperature T, so that 
p = (?p with constant sound speed c. By scaling, we assume without loss of generality that 
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with 
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as 


p = p, i. e., the speed of sound is 1 . 

Ut + F{U), + G{U)y = 0 
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Using (|^, we also write it as a symmetric hyperbolic system^ 

Ti{U)Ut + HU)Ux + B(U)Uy = 0 

with U = {p,Vi,V 2 MMy, D(U) = diag(l/p,p,p, 1 , 1 ), and 
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Applying the chain rule, we rewrite ([^ 

Ut + A{U)U, + B{U)Uy = 0 , 

where 


A = TD-^AT h B = TD ^BT 


with 


Note that, as we have used (|^ on the way from Q to (|^, the matrices A and B in (|^ 
not the Jacobians of the fluxes F and G. 



( 6 ) 

(7) 

( 8 ) 
are 


^This issue has an interesting history, cf. ElE]- 
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2 Slow and fast, parallel and non-parallel shock waves 


Ideal MHD shock waves, in their prototypical form, have the structure 
U{t,x,y) = 


U = {p ,p V ,H ), {x,y)-N<st, 

t/+ = (p+,p+V+,hf+), {x,y)-N>st, 


where N = {Ni, N 2 ) G is the direction of propagation and s the speed of the shock wave. 
Function (|^ being a weak solution of Q is equivalent to the Rankine-Hugoniot conditions 

-s{U+ - U-) + A^i(F(f/+) - F{U-)) + N 2 {G{U+) - G{U-)) = 0. 

Due to rotational and Galilean invariance it is without loss of generality that we henceforth 
assume that 

N = (1, 0) and s = 0; 

i. e. we exclusively consider shock waves of the form 


U{t,x,y) = 


U = {p ,p V ,H ), X < 0, 
t/+= (p+,p+l/+,hf+), x>0. 


and the Rankine-Hugoniot conditions read 

F{U-) = F{U+). 


(9) 


( 10 ) 


Note now hrst that for waves as for any solutions of Q whose spatial dependence is 
only via x, the divergence-free condition ([^ reduces to 


bi = a, a any constant. 

We assume 0 and simply write b,v,w instead of & 2 W 1 W 2 - 

In this paper, we are interested in Lax shocks. Following illllElIl] , two states 
U~ = {p~, p~v~, p~w~,a,b~) and U~^ = {p'^, p'^v~^, p'^w~^,a,b'^) 


( 11 ) 


that satisfy the Rankine-Hugoniot conditions (10) constitute a 

slow Lax shock iff 0 < p'^{v^)‘^ < p~{v~)‘^ < 

and a 


fast Lax shock iff < p{v~^Y < p {v )^. 


( 12 ) 


(13) 


(14) 


Two states (12) do satisfy the Rankine-Hugoniot conditions (10) if and only if the two 
quadruples {p~,v~,w~,b~) and (p"*", n"*", 6+) have coinciding images under the mapping 

^ ' ■ ■ ■ (15) 


fp\ 

V 

w 

\bj 


( pv \ 

pv'^ + p+ Ifc^ 
pvw — ab 
\ bv — aw y 


3 










that F induces by omitting its forth, trivial component, in other words if both quadruples 
satisfy the four equations 

pv = m (16) 

pv‘^ + p+ ^6^ = j (17) 

vb — aw = c (18) 

mw — ab = d (19) 

for the same values of the four parameters m,j,c,d G M. As simple argument^ show, we 
lose no generality in assuming that 

m > 0, d = 0, and pv"^ ^ a^. (20) 

Using (19) in (18) and inserting the result and (16) in ( |I^ then yields 

me 


nmc! N l+V^ 1 

g {v) = m -h - 


2 V mv — 


= j, to be solved for u G (0, a"^/m) U (a"^/m, cx)). (21) 


As for every solution u of (21), relations (16), (18), (19) provide unique associated values for 
v^w and b, understanding (21) will give a complete picture. One distinguishes two cases. 


c = 0.' parallel shocks. In this case, (21) has two solutions 


^ f 3 




J \2 


'm 


if — > 2, with 0 < < u . 

m 


The corresponding states (12) constitute a 

slow parallel shock iff mv~ < 

and a 

fast parallel shock iff < mv^. 


( 22 ) 


(23) 


The fact that the value of a has no influence on the p, u, w, b components of parallel shocks is 
easily understood by noticing that they have b = w = 0 and thus are purely gas dynamical. 


c 7 ^ 0.- non-parallel shocks. In this case, tends to 00 not only for u \ 0 und v 00 , 
but also for v —)■ c? jm. Thus for every 

3 > m, c) = min 


(21) has two solutions 

u+(a,m,c, j) < v~{a,m,c,i) < a^/m 
that consitute a slow shock. Similarly, for every 


there are two solutions 


that dehne a fast shock. 


3 > 3Lni(^,m,c) = min 

(a^/m,oo) 


a^/m < Uj (a, m, c,j) < Vj (a, m, c,j) 


^reversing, shifting, scaling, and the observation that cases with m = 0 or pv^ = a? give no Lax shocks 


4 



















3 Lopatinski determinant and critical manifold 


According to Majda’s theorjj^ [Till ITT] on the persistence of shock fronts, the local-in-time 
stability of the planar discontinous wave (|^ is determined by the behaviour of the Lopatinski 
determinant 


A : 5+ = (C+ X M) \ {(0,0)} ^ C, A(A,a;) := det(i?-(A, w), J(A, w), i?+(A, a;)), (24) 

where C+ = (A G C : JIA > 0}. While uniform stability corresponds to the non-vanishing of 
A on all of S'+, shocks with 

0 ^ A-^(O) C * X M or 0 ^ A-^(O) n (C+ x M) (25) 

are neutrally stable or strongly unstable, respectively. The ingredients of the Lopatinski 
determinant are 

J{X,uj) := X{U~^ — U~) + iuiGiU^) — G{U~)), “jump vector”, 

R~{X,oj), base of the stable space of := (AI -|- iC(;B^)(A~)“^, 

R^{X,uj), base of the unstable space of L+ := (AI -f ia;B+)(A^)“^, 

where A^, B^ denote A(t/^), B(I7^). The theory of hyperbolic initial-boundary value prob¬ 
lems mm implies that are well-dehned bundles of constant dimension. To be precise, 
it is on C+ X M = 5*+ \ (iM x M) that the Lopatinski matrices have constantly trivial 
neutral spaces and thus “consistent splitting”, i. e., stable and unstable spaces of constant 
dimensions, so that in particular 

d~ = dim span R~{X,u) and R = dim span R~^{X,u) 

are constant; for points (A,a;) G S'+ with purely imaginary values of A, the i?^(A,a;) are 
dehned as limits from the interior of 5*+ [7]. From the one-dimensional ‘Lax counting’ of 
characteristic speeds Bill. we know that 

d~ = 1 and d~^ = 3 for slow MHD shocks, (26) 

while 

d~ = 0 and d'^ = A for fast MHD shocks. (27) 

The Lopatinski determinant A being degree-one homogeneous in (A,a;), we from now on £x 
the transverse wave number to 

u = ± 1 . 

To avoid abundant notation, we also £x from now, again without loss of generality. 


p- = 1 (28) 

and use the two parameters p+, c instead of the three paramters j,m, c. For parallel shocks, 
our choice (28) implies 

v~ = \/(R, = 


^This is what our passing, in Sec. 1, through the symmetric hyperbolic formulation ([^ is needed for. 
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In this paper we concentrate on slow shocks]^ The following is a key observation. 
Theorem 1. For slow parallel MHD shocks in 0>ii with ([^ and p = 1, 

A(0,±1) = 0 zf = (29) 

+ 1 


Proof. Interesting manipulations show that one can take 



a(p+-l) 
2 0 

0 0 

0 2 

\ 0 0 


0 \ 
0 

-ay/p+ 

0 

1 / 


Together with J = (0, 0, i, 0, 0), this yields 

A(0,1) = 2i[p+(a2 + l)-(a2 + 2)]. 


□ 



Figure 1: Slow parallel shock with p 
condition (22): p~ -{v-f <a^. 


1, cf. Theorem 1. The black boundary is the Lax 


4 Symmetry breaking 


The situation of parallel shocks is degenerate as it possesses a reflectional symmetry in 
the transverse (?/-)direction. For the Lopatinski determinant this symmetry means that 

^Cf. Trakhinin’s paper [12] (and also [21 [5]) for other results. 
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A(A, —u) vanishes exactly if A(A,a;) does. Perturbing the parameter c away from 0 breaks 
this symmetry, and the zero of A that we found, for c = 0 at A = 0, splits. 


For all values of p'^, a, c that permit a (then unique) slow MHD shock wave, we write 

for the corresponding Lopatinski determinant. Starting from Theorem 1, we found the 

following. 


Numerical Observation 1. There are an e > 0 and two functions, 


R{a, c), even in c and with R{a, 0) = 


a^ + 2 
a2 + l 


, and ■j{a,c), odd in c, 


both defined on fie = {{cl, c) : a > amin{c), —e<c<e}, such that 

A^(“’‘^)’“’'^(±i 7 (a, c), ±1) = 0 for all c G (—e, e). 





Figure 2: Curves a h-)■ (i?(a, c), 7 ( 0 , c)) for some values of c between —0.01 and 0.01. The 
red curve corresponds to c = 0 and thus to the red curve in Fig. 1. 

A detailed description of the numerics is postponed to a later publication. 
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5 Emergence of unstable modes 


Do unstable modes emerge in families of shock waves that correspond to parameter val¬ 
ues which cross the critical manifold? The following is what we conclude from numerical 
computations. 


Numerical Observation 2. There are a6 > 0 and a smooth function a+ift : [0, 5) —)■ C 

with 


and 


a(ao,c, 0)=0 and /3(ao, c, 0 ) = 7 ( 00 , c) 


a(ao,c,0>0 forf>0 


such that 


A^(“0’'=)’“o+b'=(«(ao,c,0 ±^/9(ao,c,0,±l) = 0 for all f e [0,5). 


This means that for ^ > 0, 
is an unstable eigenvalue for oj 


X = a ± i fi 

± 1 . 



Figure 3: Curves ^ (oq -l- R{ao, 0), a(ao, 0, .^)) for various values of Oq. 


A detailed description of the numerics is again postponed to a later publication. 





































Remark. Both from a physics perspective and as the Evans function for non-ideal shock 
waves is intimately related to the Lopatinski determinant for their non-ideal counterparts 
iia, one expects the gallopping instability described in this paper to occur also in the presence 
of viscosity and and electrical resistivity. 
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